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We propose a method to estimate the bias voltage dependent magnetic contrast observed on 
constant current spin-polarized scanning tunneling microscopy (SP-STM) images from single point 
differential conductance (dl/dV) tunneling spectra. Depending on the number of single point mea- 
surements, the bias positions of magnetic contrast reversals and the maximally achievable magnetic 
contrast can be determined. We validate this proposal by simulating SP-STM images on a com- 
plex magnetic surface by employing a recently developed approach based on atomic superposition. 
Furthermore, we show evidence that the tip electronic structure has a major effect on the mag- 
netic contrast. Our theoretical prediction should inspire experimentalists to considerably reduce 
measurement efforts for determining the bias dependent magnetic contrast on magnetic surfaces. 
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Achieving ultrahigh information density in a controlled 
way on surfaces of materials [H HI is one of the ultimate 
goals of magnetic research nowadays for the purpose of 
future data storage technological applications. This can 
be established by the reduction of the size of magnetic 
information storage units going down to the nanoscale or 
even to single atoms Reading and writing informa- 
tion routinely from and to such magnetic units is a great 
challenge. Spin-polarized scanning tunneling microscopy 
(SP-STM) [H, employing a magnetic tip, proved to be 
extremely successful for studying magnetism on surfaces 
in high spatial resolution. Recent experimental advances 
using this technique allow the investigation of complex 
magnetic structures (frustrated antiferromagnets, spin 
spirals, skyrmion lattices, etc.) [!, 0-0] • 

In the most routinely used constant current mode of 
SP-STM, the apparent height difference between differ- 
ently magnetized surface atoms allows the discrimina- 
tion of the individual atomic magnetic properties and 
the mapping of the magnetic structure. This apparent 
height difference is called magnetic contrast. Finding 
the maximal magnetic contrast for a given surface-tip 
combination is crucial for a more efficient magnetic map- 
ping. This can be done by using magnetic tips with large 
spin polarization, or, by choosing the appropriate bias 
voltage. A very few works focused on the investigation 
of the bias dependent magnetic contrast so far. Among 
those, a magnetic contrast reversal has been reported in 
two different magnetic systems [1, 0] , and the effect was 
related to the surface electronic structure rather than to 
the effect of the tip. In another work, such contrast re- 
versals were observed during the scanning with the STM 
tip at fixed bias, and this was identified to be due to 
the magnetic switching of the tip [lOj . Moreover, un- 
der certain circumstances, a giant magnetic contrast can 



be obtained, and this effect was explained by chemically 
modified STM tips 

Distinguishing atoms with different magnetic proper- 
ties on a complex magnetic surface can successfully be 
performed by spin-polarized scanning tunneling spec- 
troscopy (SP-STS) as well. For example this technique 
has recently been used to read out information from 
an atomic scale all-spin- based logic device [l2|. Here, 
through scanning the surface with a magnetic tip, the 
measured differential conductance (dl/dV) values vary 
depending on the magnetic properties of the underly- 
ing surface atom. For such a spectroscopic detection of 
atomic magnetism, the contrast, that is the dl/dV dif- 
ference above individual atoms, should also be tunable 
by changing the bias voltage [131 ] - 

In this work, we propose a method to estimate the bias 
dependent magnetic contrast observed in constant cur- 
rent SP-STM experiments from single point differential 
conductance tunneling spectra. The key quantity for this 
estimation is the magnetic contribution to the differential 
conductance, dlMAGN / dV , that can be obtained from 
measuring tunneling spectra with oppositely magnetized 
tips. Let us assume that one can measure differential 
conductances dip /dV(z\, V) and dI^ P /dV{z\,V) at the 
same tip position above surface atom A at z\ tip-sample 
distance, where P and AP denote parallel and antipar- 
allcl tip magnetization directions, respectively, compared 
to a predefined direction. From these two measured spec- 
tra, the magnetic contribution to the differential conduc- 
tance can be calculated as 
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Using this quantity, the magnetic current can be deter- 
mined at arbitrary V bias voltages via integration: 



'Electronic address: palotas@phy.bme.hu 



1 MAGN 



(zx,V) 



d y dI MAGN ( yy 
dV 



(2) 



2 



Here, we assumed zero temperature. For a given surface- 
tip combination, by varying the upper integral limit V, 
we can determine the bias voltage, where the magnetic 
current is zero. This corresponds to a magnetic contrast 
reversal on the constant current SP-STM image, as will 
be demonstrated later. This bias position can be found 
using magnetic tips having arbitrary magnetization di- 
rection, however, the contrast can further be enhanced 
by setting the ti p m agnetization direction properly as 
discussed in Ref. Thus, from P and AP single point 
spectra measured at tip position (A, z{), the bias volt- 
age(s) can be determined, where a magnetic contrast re- 
versal occurs. 

Similarly to the magnetic current, the topographic cur- 
rent can be calculated from the same single point P and 
AP spectra as 
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and the total current with P tip alignment is 
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The same procedure can be repeated above atom B, i.e. 
P and AP single point spectra have to be measured, and 
the corresponding currents can be calculated similarly to 
atom A, e.g. I B {z\,V) denotes the total current with 
P tip alignment at the tip-sample distance Z\ and bias 
voltage V. 

In an SP-STM image of a complex magnetic surface 
the magnetic modulation is superimposed on top of the 
topographic image. Generalizing the relation between 
constant current and constant height STM images re- 
ported in Ref. [l4| to SP-STM images, we can define the 
contrast between atoms A and B on the surface at the 
tip-sample distance z\ and bias voltage V as 
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where the average tip-sample distance has been gener- 
alized to the topographic image, i.e. the constant cur- 
rent contour of the nonmagnetic part of the tunneling 
current, Itopo- We also took advantage of the expo- 
nentially decaying character of the tunneling current, 
Itopo(z,V) = I {V)e-' 2Ko ^ z . The contrast in Eq.@ 
corresponds to the apparent height difference of atom A 
and B on a proper constant current contour, and can be 
decomposed to topographic and magnetic parts, 
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If the topographic contrast Az^q PO is zero, then the to- 
tal contrast in Eq.((6]) corresponds to the magnetic con- 
trast. This happens, e.g., on surfaces consisting of chem- 
ically equivalent but magnetically inequivalent atoms. In 



this case Itopo = Itopo 



Itopo- ^ n or der to be able 



to calculate the so far unknown decay Ko(V), a third pair 
of single point P and AP spectra is needed at a different 
tip-sample distance z 2 , e.g., above atom A. Since Itopo 
can be calculated from the spectra according to Eq.Q 
at both tip-sample distances, and the tip displacement 
z 2 — z\ can be measured in experiments, the decay can 
be determined as 

^ V) = 2jz^) ' (8) 

Thus, by measuring P and AP tunneling spectra in each 
of the three well-defined tip positions above the surface: 
(A, Z\), {A, z 2 ), and (B, zi), the bias dependent magnetic 
contrast between atom A and B at the tip-sample dis- 
tance Zx, Azab(zi, V), can be predicted following Eq.fjH). 
From this function the bias position can be identified, 
where the maximal magnetic contrast can be achieved. 
Since the magnetic contrast is assumed to be exponen- 
tially decaying with respect to the tip-sample distance 
similarly to the corrugation amplitude of nonmagnetic 
STM images, the identified bias position of the maximal 
magnetic contrast should not change by recording SP- 
STM images at different constant current conditions, i.e. 
at different tip-sample distances. 

By additionally recording P and AP tunneling spectra 
at the tip position (B 7 z 2 ), Az AB {z 2 ,V) can also be cal- 
culated. Thus, measured spectra at four well-defined tip 
positions above the surface: (A, z\), (A, z 2 ), (B,zi), and 
(B,z 2 ) are, in principle, sufficient for determining the 
bias dependent magnetic contrast between atom A and 
B at an arbitrary tip-sample distance z in the tunneling 
regime: 



Az AB {z, V) = [Az A b(z!, V)] =2-1 [Az AB (z 2 , V)} -2-1 . 

(9) 

Note that since the absolute tip-sample distances are un- 
known in experiments, the magnetic contrast can be de- 
termined at distance z relative to z\ and z 2 only. 

In order to demonstrate the applicability of the pro- 
posed procedure for estimating the bias dependent mag- 
netic contrast from single point tunneling spectroscopy 
data, we perform simulations on a sample surface with 
noncollinear magnetic order. One monolayer (ML) Cr 
on Ag(lll) is a p rototype of frustrated hexagonal anti- 
ferromagnets [lj]- Due to the geometrical frustration of 
the antiferromagnetic exchange interactions between the 
Cr spin moments, its magnetic ground state is a non- 
collinear 120° Neel state [15(. Taking the spin-orbit cou- 
pling into account, two Neel states with opposite chiral- 
ities can form, and one of them is energetically favored 
1- 

We performed fully noncollinear electronic structure 
calculations on the Cr/Ag(lll) system, based on the den- 
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sity functional theory (DFT) within the generalized gra- 
dient approximation (GGA) and using the projector aug- 
mented wave (PAW) method, implemented in the Vienna 
ab initio simulation package (VASP) PjH2CJ. The com- 
putational details are reported elsewhere [9[. We used 
two tip models: one is an electronically flat, maximally 
spin-polarized (Pt = +1) ideal magnetic tip, and the 
second is a blunt Ni tip, i.e. a Ni adatom placed above 
the hollow position of a Ni(llO) surface. More details on 
the chosen tip models can be found in Ref. [l3| . 

Taking the calculated electronic structures, the dif- 
ferential conductances can be obtained following Ref. 
plj . On the other hand, the magnetic current can di- 
rectly be determined from the sample and the tip elec- 
tronic structures at the tip apex position Rtip above 
the sample surface using the atom superposition method 
0, EH Ull HI1 based on the spin-polarized Tersoff- 
Hamann model [HI, HU, [25[ , assuming elastic tunneling: 
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Here, the sum over a includes a sufficient number of sur- 
face atoms [HI with position R a and projected magnetic 
density of states (magnetic DOS), mfg. Similarly, m T de- 
notes the magnetic density of states projected onto the 
tip apex atom. Ef and Ej, are the Fermi energies of the 
surface and the tip, respectively. Moreover, the expo- 
nential expression is the tunneling transmission assum- 
ing spherical decay of the electron wave functions and an 
effective rectangular tunnel barrier, where <ps and 4>t are 
the surface and tip electron work functions, respectively. 
The shape of the electron orbitals involved in the tunnel- 
ing [2(| are neglected (independent orbital approximation 

Figure [1] shows the simulated magnetic current, 
Imagn{V), according to Eq. (fT0| z = 3.5 A above each 
Cr atom (Crl, Cr2, Cr3) in the (V3 x magnetic unit 
cell shown in the left part of Figure measured with the 
two considered tip models. For the negative bias range, 
the integral limits are reversed in order to obtain positive 
total current values, and the magnetic current has been 
calculated accordingly. The tip magnetization direction 
is fixed parallel to the Crl magnetic moment. Due to 
the frustrated magnetic structure, the Cr2 and Cr3 mag- 
netic moments have an angle of 120° with respect to the 
tip magnetization direction. Following this, the magnetic 
currents measured above the Cr2 and Cr3 surface atoms 
are equal and their values are cos(120°) = —0.5 times 
the magnetic current above Crl. Moreover, it is clearly 
demonstrated in Figure Q] that the magnetic properties 
of the tip affect the magnetic current considerably. 

The magnetic contrast is reversed at Imagn{V) = 0. 
In order to find out the corresponding bias positions for 
the two considered tips, we have to zoom in the region [0 



^ 1 ' 














s 














- \ 






/ Ml 

/ N 1 
/ l l 








/ d 
/ * 








7 


4 O^i 




— "■ " 




\ 3 


— 




^^^^ y 








- 

-> 












\ * f 




— Crl, Ideal tip 










- Cr2, Cr3, Ideal tip 
Crl. Ni tip 






1 , 1 




-- Cr2, Cr3, Ni tip 






-2 -1 




1 




2 




Bias Voltage [V] 







FIG. 1: (Color online) Simulated magnetic current, 
/iUAGiv(V'), according to Eq. lfTU)) z = 3.5 A above each Cr 
atom (Crl, Cr2, Cr3) in the (^3 x x/3) magnetic unit cell 
(see left part of Figure [3| , measured with an ideal electroni- 
cally flat maximally spin-polarized tip (gray), or a model Ni 
tip (black). The tip magnetization direction is fixed parallel 
to the Crl magnetic moment in Figure^] The inset shows the 
region between V and 1 V zoomed in. The sign change of 
Imagn(V) occurs at 0.94 V and 0.74 V for the ideal and the 
Ni tip, respectively. These correspond to magnetic contrast 
reversals. 



V, IV]. The indicated rectangular area is shown in the 
inset of Figure [T] It is clearly seen that the sign change 
of Imagn{V) occurs at 0.94 V and 0.74 V for the ideal 
and the Ni tip, respectively. Although the reversal is, in 
principle, due to the surface electronic structure [1, Q, 
the tip plays a crucial role as well, since its electronic 
structure modifies the bias position of the reversal. Note 
that the determination of the above bias positions does 
not depend on whether the tip is placed above Crl, Cr2, 
or Cr3 atoms. Moreover, we find that the bias voltage 
of the contrast reversals is stable within ±0.01 V placing 
the tip above other surface positions. 

In order to estimate the bias dependent magnetic con- 
trast, Eq.© has to be employed. The topographic 
current Itopo can be calculated similarly to Eq. (fT0]) 
by replacing the magnetic DOS by the corresponding 
charge DOS of the sample and the tip @. The decay 
of Itopo{ z ,V) can be calculated using Eq.©. In the 
simulations we determined kq(V) from Itopo( z ,V) by 
linear regression taking 150 z values in the range [0.01 A, 
7.95 A] at each chosen bias voltage with 0.01 V steps in 
the [-2.5 V, 2.5 V] interval. Figure [5] shows the bias de- 
pendent magnetic contrast between Crl and Cr3 atoms 
at two different tip-sample distances using both tip mod- 
els. The sign of the curves corresponds to the ones de- 
noted by solid line in Figure [U and the contrast reversal 
is obtained at Azc r i-Cr3(z, V) = 0. However, due to the 
denominator of Eq.©, we obtain qualitatively different 
curves compared to Figure [T] From the calculated func- 
tions, bias voltages can be identified in both the negative 
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FIG. 2: (Color online) The bias dependent magnetic contrast 
between Crl and Cr3 atoms Azc r i-Cr3(z, V) at two differ- 
ent tip-sample distances, z — 3.5 A (solid lines) and z = 4.5 
A (dash-dotted lines) calculated using Eq.© measured with 
the ideal magnetic tip (gray) and the Ni tip (black). The lo- 
cal absolute maxima of the magnetic contrasts and their bias 
values are explicitly shown in both the negative and the pos- 
itive bias ranges. Vertical dashed lines denote the contrast 
reversals, see also Figure [JJ For comparison, color dashed 
lines show the apparent height difference between Crl and 
Cr3 atoms calculated from constant current SP-STM images, 
see text for details. 



and positive bias ranges, where a local absolute maxi- 
mum magnetic contrast can be achieved. As can be seen, 
these highly depend on the properties of the magnetic 
tip. In the studied bias interval local absolute maximum 
contrasts are expected at -1.40 V and 2.50 V for the mea- 
surement with the ideal tip, whereas these bias values are 
considerably modified to -0.84 V and 1.31 V by using the 
Ni tip. The maximum absolute contrasts are obtained 
at the corresponding negative bias values for both tips. 
Moreover, as can be seen, the change of the tip-sample 
distance does not modify the bias positions of the abso- 
lute contrast maxima. We tested it for other distances as 
well, and the bias positions were found within ±0.02 V. 

For comparison, color dashed lines in Figurc[5]show the 
apparent height difference between Crl and Cr3 atoms 
calculated from constant current SP-STM images. The 
current contour was chosen in such a way that the ap- 
parent height of the Crl atom is 3.5 A at all considered 
bias voltages. The agreement with the estimated mag- 
netic contrast using Eq.(|6]) is remarkable, particularly for 
Azcri~Cr3 > 0. There are some quantitative deviations 
if Azcvi-cv3 is negative, but the overall qualitative bias 
dependence is well reproduced by Eq.©. 

In order to better visualize the bias dependent con- 
trast, we simulated constant current SP-STM images. 
Figure [3] shows such images measured with the ideal mag- 
netic tip (top row) and the Ni tip (bottom row) at V, at 
the bias voltages corresponding to the contrast reversal 
(0.94 V and 0.74 V, respectively), and at the voltages cor- 



FIG. 3: (Color online) Simulated SP-STM images depending 
on the bias voltage and the considered tip: ideal magnetic tip 
(top row), Ni tip (bottom row). The magnetic moment of the 
tip apex atom is parallel to Crl, and indicated by a vector 
(M TIP ). The surface geometry of 1 ML Cr on Ag(lll) and 
its ground state magnetic structure are explicitly shown. The 
Cr and Ag atoms are denoted by spheres colored by green 
(medium gray) and purple (dark gray), respectively, while 
the magnetic moments of individual Cr atoms are indicated 
by (red) arrows. The {y/3 x y/3) magnetic unit cell is drawn 
by yellow (light gray) color, and the scanning area is shown 
by the black-framed rectangle. 

responding to the local absolute maxima of the magnetic 
contrast in both the negative and positive ranges. The 
surface geometry is shown in the left part of the figure, to- 
gether with the scanning area (black- framed rectangle), 
and the magnetic unit cell (yellow (light gray) framed 
rhombus). The tip magnetization direction is parallel to 
the Crl magnetic moment. For the zero bias images a 
temperature of 4.2 K has been considered following Rcf. 

0, since there is no current at K. 

We find a similar type of magnetic contrast at V and 
at the negative bias voltages for both tips, respectively, 

1. e., for V < 0, Azc r i-Cr3 > (Crl appears higher than 
Cr2 and Cr3) for the ideal tip, and Azcvi-cv3 < (Crl 
appears lower than Cr2 and Cr3) for the Ni tip, similarly 
as reported in Figure [5] At the corresponding reversal 
voltages, 0.94 V (ideal tip) and 0.74 V (Ni tip), all Cr 
atoms appear to be of equal height on the SP-STM im- 
age. Here, the magnetic contrast is lost as Imagn = 0. 
This corresponds to a qualitatively similar image of per- 
forming the STM measurement with a nonmagnetic tip 
@. Thus, the magnetic current calculated at a single 
point above the surface was indeed able to find the cor- 
rect bias position of the magnetic contrast reversal. At 
biases above the reversal voltage, the magnetic contrast 
is inverted. We illustrate this for the image calculated at 
2.50 V for the ideal tip, and at 1.31 V for the Ni tip. 

In conclusion, we proposed a procedure for determining 
the bias dependent magnetic contrast from single point 
differential conductance tunneling spectra measured with 
oppositely magnetized tips. By performing simulations 
based on first principles electronic structure calculations, 
we demonstrated that the obtained bias dependent con- 
trast corresponds well to the apparent height difference 
calculated from constant current SP-STM images. Our 
theoretical prediction should inspire experimentalists to 
considerably reduce measurement efforts for determining 
the bias dependent magnetic contrast on magnetic sur- 
faces. 
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